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Abstract 



"^ ■ The mathematical notion of a spectral singularity admits a physical interpretation as a 

"i ' zero-width resonance. It finds an optical realization as a certain type of lasing effect that 

"jrt I occurs at the threshold gain. We explore spectral singularities of a complex spherical barrier 

potential and study their realization as transverse spherical electromagnetic waves emitted by 

a gain medium with a spherical geometry. In particular, for a typical dye laser material, we 

^ ' obtain a lower bound on the size of the gain medium for the occurence of this kind of spectral 

singularities. 



PACS numbers: 03.65.Nk, 42.25.Bs, 42.60.Da, 24.30. Gd 
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An interesting by-product of the recent study of complex potentials supporting a real spectrum 
[H 12] is the discovery of the optical realizations [21 HI El E] of the mathematical notion of a spectral 
singularity [7]. For a complex scattering potential defined on the real line, these are certain points 
of the continuous spectrum of the corresponding non-Hermitian Schrodinger operator at which the 
reflection and transmission coefficients diverge [3]. As a result, they correspond to scattering states 
that behave exactly like resonances. This observation has also found applications in condensed 
matter systems [5] and triggered further study of the subject ^. 

The study of the optical realizations and applications of spectral singularities has so far been 
confined to effectively one-dimensional models involving infinitely long waveguides [3l H] or infinite 
planar slab gain media [HE]- In the present article we examine spectral singularities of a three- 
dimensional complex spherical barrier potential that admits a physical realization involving an 
optical gain medium with a spherical geometry. 
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Consider the propagation of electromagnetic (EM) fields in a linear source-free isotropic dielec- 
tric with time- independent (complex) permittivity 5 = e(r^ and constant permeability fi = fiQ. The 
electric and magnetic fields, E and B, satisfy Maxwell's equations 

V-Z3 = 0, V-5 = 0, (1) 

V X E + B = 0, V X H -3 = 0, (2) 

where D := e{r)E, H := fioB, and a dot represents the time-derivative. Following the standard 
derivation of the wave equation, we can use ([2]) to obtain 

E + n'^E = (3) 

where Q"^ := [yUoe(r)]^^V x V. 

For a time-harmonic EM field propagating in a spherically symmetric dielectric we have E{f, t) = 
e~^^^E{r,t) and e{r) = Bo^{r), where r is the radial spherical coordinate and 
3 : [0, oo) — )• C is a complex-valued function. In this case ([3]) reduces to the time-independent 
Schrodinger equation: 

- V^E{r} + d{r)E{r) = k'^E{r), (4) 

where k := w/c and d{r) := /c^[l — 3(r)]. The latter is a spherical complex barrier potential provided 
that the dielectric medium is confined in a spherical region. This is the case for 

f n^ for r < a, f P(l - n^) for r < a, 

I 1 tor r > a^ I U tor r > a, 

where n is the complex refractive index of the dielectric medium that is supposed to be constant. 
Consider spherically symmetric solutions of (j4]) that define transverse radially propagating 
spherical waves. They have the general form E{f) = Eg{r)9 + E^{r)(j), where 6 and are the 
unit vectors along the 6- and 0-directions. Inserting this ansatz in (^ we find that either 6 = 7r/2 
or Eq = 0. In what follows we consider the latter case, so that E{r) = E^{r)(f), and (jl]) reduces to 

-^ + -— + /;:'- ^(r)-- 
dr^ r ar r^ 

For r < a and r > a where 'd takes constant values, we can transform ([6]) to the Bessel equation 

of order u := v5/2, [inj. The electric field inside and outside the dielectric region are respectively 

given by 

^outside = [A,h<;}\kr) + A2h<;^\kr)]^, (7) 

inside = [Bij^ikr) + 52n^(^r)]0, (8) 

where the coefficients Ai^2 and Bi^2 are related via the appropriate boundary conditions at r = a, 

(1 2) ~ 

jy and hu ' are respectively the spherical Bessel and Hankel functions, and k := n k. Demanding 
the electric field to disappear at the origin, we have B2 = and 

Einside = Bijy{kr) 0. (9) 



E^ = 0. (6) 



In view of the asymptotic expression for the Hankel functions, namely 
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the terms on the right-hand side of ([7]) correspond to the reflected and incident spherical waves, 
respectively. As a result, we identify R := |yli/y42p with the reflection coefficient of the system 
whose real and positive poles in the /c^-plane are the spectral singularities. In order to determine 
these poles, we need to derive the relationship between the coefficients A12 and Bi. First, we use 
(j2]) to compute the magnetic field: 
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Imposing the standard boundary conditions: (-E'||)insidc 
noting that H = fioB, we have 



(£'11 ) outside) ( -f^l I ) inside — (-f^||) outside and 
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We can easily eliminate Bi in flT^ and flT^ . and compute the reflection amplitude: 
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(14) 



Setting the denominator of this quantity equal to zero we flnd the following condition for the 
existence of a spectral singularity. 
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Recalling the recursion formula [10] for hy and jy-. 

duyjz) _ vuy^iiz) - {u + 1)m^+i(z) 
dz ~ 2z/ + 1 

with Uy = jy or hy\ we can express ( IT5l) in the form 

[i^/i2i(/ca) - (z/ + l)h^^l^{ka)\ n [ujy_i{ka) - (z/ + l)jy+i{ka)] 



(16) 
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where we have used the fact that k = nk. Because n and consequently k take complex values, 
this is a complex equation. Equating the real and imaginary parts of both sides of this equation we 



find two real equations for tlie tliree real variables t] := Re( n), k := Im( n), and x := ka. Because 
the determination of the explicit form of these equations is not easy, we examine their asymptotic 
expressions that are valid for x ^ 1. For typical situations the radius of the dielectric ball is much 
larger than the wavelength of the wave. Therefore, for practical purposes that we will consider, the 
condition x ^ 1 holds and the asymptotic treatment of (IT7|) provides extremely reliable results. 

The asymptotic expansion of the spherical Bessel and Hankel functions, j,y and hi , that are 
valid in the large argument limit (|z| ^ 1), have the form 
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where 



A,M := 
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and r stands for the Gamma function. Substituting (TTSi) and (TT9l) in (TT7|) . neglecting the quadratic 
and higher order terms in x~^ and noting that Aq{v) = 1 and Ai := Ai{v) = " ^^ , we find 

n^-l)Ai 



tan(na; —J 
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(20) 



Next, we wish to compute the real and imaginary parts of ( 1201) and express them in terms of 
the real parameters r], k, and x. To this end, first we recall the identity: 



-1/ \ 1 , /I + zz 

tan (z) = nm -\ m , , 

2i \l — iz 



(21) 



where z is a complex variable, m is an integer, and "In" denotes the principal part of the natural 



logarithm. If we set z = —in 
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and use fl20|) and fl2T]) . we obtain 
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Note that 



In 



-1 + 



2iAi 



ln(-l)+ln 1 - 



2iAi 



nx J \ nx 

Therefore, if we only keep the first two largest powers of x in 
drops out of the calculations and we obtain 
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Next, we employ n = 77 + m to obtain 

1 /(r/ + l)2 + K2 



[j] - 1)2 + fi;2 



exp 
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where "arctan" stands for the principal part of "tan~^" with values in (— f , f )• Furthermore, 
because for realistic gain media k ^ 1 ^ r/, we neglect the quadratic and higher order terms in k. 
Substituting n = t] + in and flMl) in fl23|) . we arrive at 

XT] ^ vr I ?Ti H J, (25) 

.. . -i,n(i±i). (.«) 



Equivalently, 



-^ln(ff 
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(27) 



Because ?7 and x take positive values and x ^ 1, the integer m that we view as a mode number 
take large positive integer values. This implies that n < 0. Therefore, in order to create a spectral 
singularity the dielectric must consist of a gain medium. This is an intriguing mathematical mani- 
festation of the law of conservation of energy. The same results is obtained in [U [5] while studying 
the optical spectral singularities of a one- dimensional waveguide including a gain region and an 
infinite planar slab gain medium. 

Equation fl271) determines the approximate location of the spectral singularities in the rj-n plane 
(complex n -plane). It yields an infinite family of curves parameterized by the mode number m. 
Given that x = ak, for each point (r/^, k^) that lies on one of the spectral singularity curves. 
Equation (1281) specifies the value of /c^ for the corresponding spectral singularities. 

In practice the refractive index n depends on the properties of the gain medium and the 
wavelength A := 27T/k = 27rc/uj of the propagating EM wave. In order to realize an optical spectral 
singularity, we need to impose the dispersion relation that determines the dependence of n on A 
and other relevant physical parameters. 

Following [5] , we consider a gain medium that is obtained by doping a host medium of refraction 
index uq and modeled by a two-level atomic system with lower and upper level population densities 
Ni and N^, resonance frequency Uq, and damping coefficient 7. Then the dispersion relation takes 
the form 

n' = nl- ^^ ;^^^.. . , (29) 

where u := oj/uq, 7 := •y/ujQ, Cjp := {Ni — Nu)e^ / {171^6^00 q) , e is electron's charge, rrie is its mass, 
and Eq is the permittivity of the vacuum. We can express Cj'i in terms of the imaginary part kq of 
n at the resonance wavelength Aq := 2iic/ujq according to [5] 

ul K. 2no7«o, (30) 



where the approximation symbol means that we neglect quadratic terms in kq. 
Inserting (l30|l in (129|) and using n = 77 + m, we obtain 

r]^no + Kofi{'^,u), KRi ^0/2(7,^), (31) 

where 

r /. .X 7(1-^2) 7^0; 

(1 — w"')"' + 7 w (1 — w"^)^ + 7"^a;^ 

Now, we substitute ( l3T|) in ( l25l) . recall that x = ka = 27ra/A, and keep the leading order term in 
the small kq approximation. This yields 

A - '"°" ■ (33) 

Next, we recall that kq may be related to the gain coefficient go at the resonance frequency according 
to go = —A-kko/Xq. Using this relation and Equations ( 125|) and (13T]) . we obtain 

^°'^Ao(2m + z. + l)/2(7,Ao/A)- ^^^^ 



In view of ( !33l) this relation takes the form 



1 , /^?^o 
5-0 ~ - In 



1' 



/(A) := ( ^^ ) . (35) 



These equations describe the dependence of the gain coefficient go necessary for generating a spectral 
singularity on the radius a of the spherical gain medium, the normalized damping coefficient 7, and 
the wavelength A. go is inversely proportional to a. Its dependence on A is given by the function 
/ which has a minimum at the resonance wavelength Aq and increases monotonically as |A — Ao| 
increases. This shows that the spectral singularity with the wavelength closest to the resonance 
wavelength of the gain medium requires the smallest necessary gain. This is in complete agreement 
with the results of [B]. 

Because our model involves no mirrors, we need a high-gain medium to create a spectral sin- 
gularity. In the following we consider a semiconductor diode laser (disregarding the difficulty of 
manufacturing a spherical diode) and a dye laser. 

Consider a spherical diode laser with the following characteristics jllj : 

no = 3.4, Ao = 1500nm, 7 = 0.02, ^o<1000cm-\ (36) 

Then according to ( 1331) the spectral singularity with wavelength closest to the resonance wavelength 
Ao corresponds to the mode number m = 679, if we take a = 150 /im. Its wavelength and the 
required gain coefficient are respectively given by A = 1499.870 nm and go = 40.412 cm"^. There 
are a total of 66 spectral singularities with mode numbers in the range 647-712 whose realization 
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Figure 1: Spectral singularities of a spherical dye gain medium with specifics ( !37|) and radius 
a = 3.300 mm. The mode number m ranges over 17746-17812 as the wavelength A decreases. The 
spectral singularity requiring the least gain coefficient Qq has a wavelength that is closest to the 
resonance wavelength of the medium. This corresponds to m = 17779. The grey horizontal line 
represents the experimental upper limit on g^. 

requires a gain coefficient less that 1000 cm~^. Their wavelength decreases from 1573.930 nm to 
1430.457 nm as m takes values from 647 to 712. 

A more realistic choice for the gain medium that supports the above-examined spectral singu- 
larities is a spherical dye laser. For definiteness consider a Rose Bengal-DMSO (Dimethyl sulfoxide) 
solution with the following characteristics [HI [12] . 

no = 1.479, Ao = 549nm, 7 = 0.062, go<5cm~\ (37) 

Using these figures in (!33l) and (134|) we find out that the smallest value of radius a that supports 
a spectral singularity (and highest gain coefficient go = 5cm~^) is 3.287825 mm which is within 
the experimentally attainable range [11]. For a slightly larger sample with radius a = 3.300 mm, 
we can generate 67 different spectral singularities by varying the gain coefficients in the range 
4.981546 — 4.999727 cm~^. These correspond to the mode numbers 17746-17812. As m takes 
values in this range, the wavelength of the corresponding spectral singularities decreases from 
550.028673 nm to 547.991700 nm. Figure [1] shows the location of these spectral singularities in 
the X-go plane. The dependence of the quantities m, A and go on one another is similar to that of 
the spectral singularities of the one-dimensional system studied in [5l [6]. As one starts pumping 
the gain medium there arises no spectral singularities unless the gain coefficient reaches its first 
critical value ((^g = 4.981546 cm~^.) This is the gain required to excite the spectral singularity 
with a wavelength (A^^-* = 549.008297507 nm) that is closest to the resonance wavelength of the 
medium. As go increases further this spectral singularity disappears and there are no other spectral 
singularities until go takes its second critical value, namely g^ = 4.981554 cm^^. This excites the 
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549.00830142 
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4.981554 
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548.97742540 


548.97743614 
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4.981572 
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17781 
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548.94657824 
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4.981630 
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549.07006387 
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548.91568378 


548.91572380 
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4.981720 


17776 


549.10095031 


549.10090243 



Table 1: The first seven critical values of the gain coefficient g^ and the wavelength of the corre- 
sponding spectral singularities of the spherical dye gain medium ( 1371) with radius a = 3.300 mm. 
Apert. ^-^d Agxict ^^6 respectively the values of the wavelength obtained using perturbation theory 
and exact numerical calculations, m is the mode number. 

spectral singularity with wavelength A'-^-* = 548.977436136 nm. This process continues until one 
saturates the upper experimental bound on the gain coefficient. As noted in |5l |6], it can be used 
to generate a tunable laser whose wavelength may be adjusted by changing the pumping intensity. 

Table [H gives the values of the first seven critical values of qq and A that we have obtained 
using perturbative and exact (numerical) calculations. The extremely good agreement between 
the perturbative and exact results is related to the fact the perturbative calculations amount to 
ignoring terms of order a;~^ and Kq. For the dye laser we consider these are about 7 x 10~^° and 
5 X 10~^°, respectively. 

Figure [2] shows a logarithmic plot of the exact expression for the reflection coefficient, R := 
|y4i/y42p, as a function of the wavelength A for the above spherical dye laser with radius a = 
3.300 mm and gain coefficient go = do ■ The highest peak with R > 5 x 10^^ demonstrates the 
spectral singularity with wavelength A = A^^^ The other peaks with i? < 4 x 10^^ are the resonances 
that turn into spectral singularities as one increases Qq. If we use more and more date points to 
produce this graph, the height of the central peak increases indefinitely while those of the others 
remain essentially unchanged. This is a clear indication that the central peak corresponds to a 
spectral singularity. 

In conclusion, we have solved the problem of finding spectral singularities of a complex spherical 
barrier potential and outlined a possible physical realization of the corresponding resonance effect. 
This is done by examining spherical transverse electromagnetic waves emitted by a spherical gain 
medium. We have obtained the particular values of the physical parameters of the system for which 
a spectral singularity appears. In particular, we have performed a reliable perturbative analysis 
of the problem that reveals the emergence of a mode number and leads to an explicit relationship 
between the necessary gain coefficient, the radius of the sample, and the wavelength of the spectral 
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Figure 2: (Color online) Logarithmic plot of the reflection coefficient as a function of the wavelength 
for the spherical dye gain medium with specifics fl37|) . radius a = 3.300 mm, and gain coefficient 



9o = % 4.981546 cm~^ that produces the spectral singularity with wavelength A = A^^'' = 

549.00329751 nm. This corresponds to the central peak. Other peaks represents resonances that 
turn into spectral singularities if we increase Qq. These have a height that is at least 100 times 
smaller than that of the spectral singularity at A = A*^^-*. The peak to the left of the central peak 
appears at A = A*-^-*. It is the first resonance that becomes a spectral singularity as go increases. 
This explains why it has the second largest height. 



singularity. For definiteness, we examined the numerical values of these quantities for a concrete 
dye laser and obtained a lower bound on the radius of the sample that would support a spectral 
singularity. This gives a ~ 3.29 mm that is much larger than the typical mirosphere dye lasers. The 
reason is that the spectral singularities we study correspond to a lasing effect that is different from 
the lasing due to the whispering gallery modes [13]. This is simply because we consider spherical 
electromagnetic waves. 
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